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^ ■ ABSTRACT 

o : 

The degree of polarization in L dwarfs of spectral types LO to L6 resulting 
c$ ! from dust scattering in a rotation-induced oblate photosphere is calculated. As- 

j> \ suming that forsterite is the main condensate, the atmospheric dust distribution 

is derived for different spectral types based on a chemical equilibrium model. 
The degree of polarization at optical is then calculated using a single scattering 
model. The expected linear polarization at optical is found to peak at around 
spectral type LI. For a fixed rotational velocity, the degree of polarization de- 
creases from hotter to cooler objects. However, with the increase in mean grain 
size, the degree of linear polarization reduces significantly. We fit the recently 
observed linear polarimetric data of L dwarfs and find that single dust scattering 
model coupled with the chemical equilibrium models of condensates is consistent 
with the observational results. 
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1. Introduction 

Observations of L dwarfs with effective temperatures of 1400-2200 K have led to the 
investigations of dust condensates in their atmospheres (Kirkpatrick et al. 1999). Be- 
cause of complete gravitational settling, grains are expected to condense beyond the visible 
atmosphere for objects with effective temperatures below 1400 K. At higher effective tem- 
peratures, grains can be present in the visible atmosphere due to incomplete gravitational 
settling (Burrows & Sharp 1999; Burrows et al. 2001; Ackerman & Marley 2001; Allard et 
al. 2001; Tsuji, Nakajima & Yanagisawa 2004; Tsuji et al. 1996). 

Using a thermodynamical model based on homogeneous and heterogeneous condensa- 
tion assumptions, Cooper et al. (2003, hereafter C03) were able to determine the distribution 
of different species of dust particles. Using solutions of dust moment equations in a static 
atmosphere, Woitke & Helling (2003, 2004) studied the continuous nucleation of solid par- 
ticles from the gas phase. By imposing the requirement of minimum mixing activities and 
thermodynamical stability, they were able to determine the upper boundary (called "cloud 
deck" ) and the lower boundary (called "cloud base" ) respectively. 

The possibility of detecting polarization at optical from grains in the atmospheres of L 
dwarfs was first raised by Sengupta & Krishan (2001). Since fast rotation of L dwarfs (Basri 
et al. 2000) will induce the shape of their photosphere into the form of an oblate ellipsoid, 
this non-sphericity will lead to incomplete cancellation of the polarization from different 
areas of the stellar surface. This prediction was first confirmed by the detection of linear 
polarization at 768 nm from a few L dwarfs by Menard, Delfosse & Monin (2002). Recently 
Zapatero Osorio et al (2004), have reported R and I band detection of linear polarization 
from several L dwarfs. Since polarization in the optical is unlikely to be due to Zeeman 
splitting of atomic or molecular lines or by synchrotron radiation, the observed polarization 
can best be explained by single dust scattering in a rotationally induced oblate atmosphere 
(Sengupta 2003). 

The observed polarization of DENIS-P J0255-4700 and 2MASSW J0036+1821 was mod- 
eled by Sengupta (2003) by assuming minimum oblateness of these objects based on a model 
of complete degeneracy of a non-relativistic polytropic gas under hydrostatic equilibrium. 
The observed degree of polarization was obtained by considering different parameters for the 
scale height, number density, and mean radius of the grains. By minimizing the oblateness 
and the number density of grains, an upper limit on the mean grain radius that is consistent 
with the polarization observed by Menard, Delfosse & Monin (2002) was obtained. However, 
the process of dust condensation was not considered. 

In the present paper, we present detailed theoretical models of optical linear polarization 
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from L dwarfs of spectral types ranging from LO to L6 based on the dust condensation model 
of C03. With this model, we fit the observed polarization from several L dwarfs and predict 
the amount of polarization expected from any L dwarfs with given rotational velocity. 

In section 2, we present the formalism for the calculation of polarization from single dust 
scattering in an oblate medium. In section 3, the adopted atmospheric model is described. 
The estimation of rotationally induced oblateness is presented in section 4, and the dust 
model parameters are described in section 5. The results and discussions are presented in 
section 6, followed by conclusions in section 7. 



If a stellar object is perfectly spherical then the net polarization would be zero due to 
the cancellation of the contribution of each point on the photosphere. The observation of 
non-zero polarization from a few L dwarfs therefore suggests that the scattering geometry 
is asymmetric, which could be the result of fast rotation of the objects. Since the dust 
density is low and scattering by atoms and molecules does not contribute to polarization 
significantly, single scattering approximation is reasonable for the region where the optical 
depth r < 1. If present, multiple scattering can reduce the degree of polarization by a few 
orders of magnitude (Sengupta & Krishan 2001) because the planes of the scattering events 
are randomly oriented and average each other's contribution out from the final polarization. 
The effects of oblateness of an object on linear polarization due to single scattering by 
spherical grains have been discussed by Dolginov, Gnedin & Silant'ev (1995) and by Simmons 
(1982). In this paper, we adopt the formalism given by Simmons (1982). 

For an optically thin atmosphere, the total linear polarization p(k) can be written as 



where k = 2n/\, A being the wavelength, the asterisk denotes the complex conjugate and 



where dw is the element of solid angle. In the above equation, 9 is the scattering angle, n 
is the number density of scattering particles, i\ and i 2 are the scattering functions given by 
(van de Hulst 1957) 



2. Polarization by single dust scattering 



p(k) = \Z(k)\ = \Z*(k)\, 
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and 



12(0) 



where 



and 
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n n (9) = -^-P^(cos9), 
sin (7 

^W = -^P„ 1 (cos^). 
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The co-efficients a n and b n are in general complex functions of the refractive index m and 
the particle radius to wavelength ratio. 



For a smooth density distribution, one can write 

00 m=l 



n(r,6,(f>) = ^2Y1 n im(r)Y lm (6,<f>), 



1=0 m =-l 

where 

Yi m {6,4>) = a(l,m)P[ n (cos 8) exp(im^) 
\ (2l + l)(l-m) ] 1/2 

and P; m is the associated Legendre function of the first kind. 
Substituting eq. 7 into eq. 2 and integrating over 0, we get 

27T ' 
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where 
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a(l, m) 



i\ ~ 12 



P{ n (cos9)d(cos9). 
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Considering an axisymmetry density distribution with a rotational invariance around 
some axis and using the addition theorem of spherical harmonic, Ni m can be written as 



N lm = 2na(l, m)P[ n (cos 9i)e 



-2i<f> 



Ri 



n(r)dr 



Pi(n)dii 



_! [1 + (A 2 -1)// 2 ]V2' 



(12) 



where R\ and R 2 are the outer and the inner equatorial axis length, A is the ratio of the 
length of the equatorial axis to the polar axis, 9{ is the viewing angle of the observer, and 

pL = COS 9. 
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At an edge-on view, Q { = n/2 and = 0, and hence Ni m is real. From the equation of 
hydrostatic equilibrium, we have an expression of the particle density as a function of gas 
pressure: 

n(P)dP , . 

" (r)dr = W' (13) 

where P is the pressure at different geometrical height, p is the mass density at different 
pressure scale, and g is the surface gravity (which can be assumed to be constant for a 
geometrically thin atmosphere). Substituting eqs. 11 and 12 into eq. 10, we have 



for the degree of polarization for a hydrostatic atmosphere viewed edge on. 



(14) 



3. The atmospheric models 

The effective temperature of the L dwarfs of different spectral type can be approximated 
by the linear empirical relationship of Stephens et al. (2001): 

T eff = 2220 — 100 x L (15) 

where L is the spectral type between L0 and L8. Recently Vrba et al. (2004) and Golimowski 
et al. (2004) have presented T e ff measurements for L and T dwarfs based on bolometric 
luminosities. In the present work we adopt the sixth order polynomial fit given by Golimowski 
et al. (2004) for translating the T e ff into spectral type. Fig. 1 shows the T e ff for different 
spectral types calculated by using the linear relationship given by Stephens et al. (2001) and 
by the polynomial formula provided by Golimowski et al. (2004). The T e // calibration of 
Vrba et al. (2004) and Golimowski et al. (2004) agree well in the interval L3-L8 but there are 
significant differences in earlier types. In our calculations for the degree of polarization, the 
effective temperature T e ff is used and hence the degree of polarization should be considered 
strictly as a function of T e ff rather than the spectral type. 

The surface gravity of L dwarfs older than about a few hundred million years varies 
from g = 10 5 cm s~ 2 to g — 3 x 10 5 cm s -2 . Evolutionary models by Chabrier et al. (2000) 
show younger L dwarfs to have surface gravity smaller than 10 5 cm s~ 2 . In this paper, we 
assume g = 10 5 cm s -2 . 
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Formation of dust makes it a prohibitive task to develop a fully consistent atmospheric 
model for ultra-cool dwarfs. This is mainly because of the fact that the presence of dust 
cloud affects the radiative equilibrium of the upper atmosphere and hence alters the T-P 
profile from that of a cloud-free atmosphere. On the other hand, the T-P profile dictates 
the position and the chemical equilibrium of condensates. Allard et al. (2001) presented 
atmospheric models for two limiting cases, e. g., one with inefficient gravitational settling 
wherein the dust is distributed according to chemical equilibrium predictions (AMES-dusty) 
and another with efficient gravitational settling in which situation dust has no effect on the 
thermal structure ( AM ES-cond) . Tsuji, Nakajima & Yanagisawa (2004) have proposed a 
Unified Cloudy Model (UCM) in which the segregation of dust from the gaseous mixture 
takes place in all the ultra-cool dwarfs and at about the same critical temperature. 

In the present work, the temperature-pressure profiles for the L dwarfs with different 
spectral type are calculated by solving the non-LTE radiative transfer equations coupled with 
the hydrostatic equilibrium equations. We have employed the full Mie theory, incorporating 
the dust opacity as well as the Mie phase function. The calculations are first performed by 
taking only gas opacities. We have adopted the atmospheric opacity sources discussed in 
Saumon et al. (2000) and assumed solar metallicity. The T-P profile thus obtained deter- 
mined the position of the cloud location. We then incorporated dust opacities to calculate 
the final set of T-P profile. This T-P profile is then used again to determine the base and 
the deck of the cloud. Our pressure-temperature profiles are checked against those presented 
in Burrows et al. (2001) and in C03, as well as models kindly provided by M. Marley (pri- 
vate communication). Nevertheless, the polarization profile would change depending on the 
chemical equilibrium procedure adopted in different atmospheric models. 



4. Rotation induced oblateness 

The observation of non-zero polarization from unresolved objects indicates the non- 
sphericity of the photosphere. The non-sphericity may results from several causes. For 
example, rotation of a stellar object will result in a shape of an oblate ellipsoid, as is evident 
in the outer solar planets. At 1 bar pressure level, the eccentricity of Jupiter, Saturn and 
Uranus are 0.35, 0.43 and 0.21 respectively. Apart from rotational effects, tidal interaction 
with the companion in a binary system also imposes an ellipsoidal shape extending toward 
the companion. 

Spectroscopic studies by Basri (1999) and Basri et al. (2000) indicate rapid rotation of 
ultra-cool dwarfs along their axis. The brown dwarf Kelu 1 is found to be the fastest rotator 
with a projected angular velocity {y shu) as high as 60 km s -1 . The observation of optical 



-7- 



polarization from L dwarfs with known projected angular velocity implies non-sphericity of 
the photosphere due to rotation. 

The oblateness of a rotating object has been discussed by Chandrasekhar (2003) in the 
context of polytropic gas configuration under hydrostatic equilibrium. For a slow rotator, the 
relationship for the oblateness / of a stable polytropic gas configuration under hydrostatic 
equilibrium is given by 

< i6 > 

where M is the total mass, R e is the equatorial radius and Q is the angular velocity of the 
object. C is a constant whose value depends on the polytropic index. For the polytropic index 
n — 0, the density is uniform and C = 1.875. This configuration is known as the Maclaurin 
spheroid. For a polytropic index of n — 1.0, C = 1.1399, which is appropriate for Jupiter 
(Hubbard 1984). For non-relativistic completely degenerate gas, n = 1.5 and C = 0.9669. 
The rotationally induced oblateness of solar planets has been discussed in details by Hubbard 
(1984) and Murray & Dermott (2000). Recently, the formalism for oblateness is extended 
to extra-solar planets by Barnes & Fortney (2003) who used Darwin-Radau relationship 
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(Murray & Dermott 2000; Barnes & Fortney 2003) to relate rotation to oblateness. In 
eq. 17, K — I /MB? e < 2/3 is the moment of inertia parameter of an object with moment 
of inertia /. The Darwin-Radau relationship is exact for uniform density objects (K = 0.4) 
and provides a reasonable (within a few percent of errors) estimation of the oblateness of 
the solar planets. 

Since L dwarfs have an extended convective region and their moment of inertia is an 
undetermined parameter, we adopt the relationship given by Chandrasekhar with the poly- 
tropic index n = 1.0 and n = 1.5. n < 1.0 would yield a too low density for brown dwarfs, 
whereas n — 1.5 would provide minimum possible oblateness due to rotation. We have cal- 
culated the mass and radius of L dwarfs of different spectral types by adopting the empirical 
relationship given in Marley et al. (1996). 



5. The dust parameters 



The dust distribution in the atmosphere is calculated based on the one dimensional 
cloud model of C03. This model assumes chemical equilibrium throughout the atmosphere, 
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and uniform density distribution across the surface of an object at each given pressure and 
temperature. Under these assumptions, the number density of cloud particles is given by 



where p is the mass density of the surrounding gas, a is the cloud particle radius, p d is the 
mass density of the dust condensates, p and p d are the mean molecular weight of atmospheric 
gas and condensates respectively. The condensate mixing number ratio (q c ) is given as 



for heterogeneously condensing clouds. In the above equation, qbeiow is the fraction of con- 
densible vapor just below the cloud base, P c j is the pressure at the condensation point 
(presented graphically in figure 1 of C03), and P is the gas pressure in the atmosphere. As 
in C03, we employ the ideal gas equation of state p = pP/RT to relate the gas mass density 
p to the atmospheric temperature T and pressure P. R is the universal gas constant and 
p = 2.36 for solar composition gas in which hydrogen is present in the molecular state. 

The kinds of solid condensates possible in the atmospheres of L dwarfs include forsterite 
(Mg 2 Si0 4 ), gehlenite (Ca 2 Al 2 Si0 7 ) and iron (C03), as well as Ti0 2 (Woitke & Helling 2004). 
However, forsterite forms in abundance in the atmosphere and plays the most crucial role 
in governing the continuum spectrum. The other species such as gehlenite are much lees 
abundant than forsterite (about a factor of 10). In the present work, we have consider only 
forsterite. The effect of other species on the degree of polarization is discussed later. We 
have used the data from C03 for forsterite : p d = 140.7 g mol -1 , p d = 3.2 g cm -3 , and 
Qbeiow = 3.2 x 10~ 5 (assuming solar abundance distribution of the elements as discussed in 
C03). Throughout our investigation, we have considered the wavelength A = 0.850 pm as 
the observed polarimetric data presented by Zapatero Osorio et al (2004) is obtained at I 
band centered at this wavelength. The object 2MASS J0036+1820 is observed by Zapatero 
Osorio et al (2004) at R Band (centered at 0.641 pm) as well as at I band (centered at 0.85 
pm). The same object is observed by Menard, Delfosse & Monin (2002) at I band (centered 
at 0.768 pm). Therefore, for this object we model the degree of polarization at different 
wavelengths. The real part of the refractive index is fixed at 1.65 and the imaginary part 
is taken by interpolating the data given in Scott & Duley (1996). It should be mentioned 
that the refractive index of amorphous condensates might differ under different physical 
conditions. 




(18) 



Qc Qbeiow 



(19) 
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Apart from the calculation of the grain number density, the location of the cloud in the 
atmosphere plays an important role in determining the amount of polarization. The location 
of the cloud base for different atmospheric models and different chemical species is determined 
by the intersection of the T-P profile of the atmosphere model and the condensation curve 
P Cy i as prescribed in C03. Taking the condensation curve for forsterite, we determine the 
base of the cloud for each spectral type, from LO to L8. Figure 2 presents the atmospheric 
pressure height at which the cloud base is situated for models with different spectral types. 
As the effective temperature decreases from LO to L8, the cloud base is pushed deeper into 
the atmosphere. According to the condensation curve given in C03, the base of forsterite 
cloud for a L8 object (T cff ~1480 K) is situated at about 10.0 bar pressure height when 
the surface gravity of the object is assumed to be 10 5 cm s -2 . For a similar model, Woitke 
& Helling (2004) found the base of TiC>2 at 9.4 bar pressure height. Similarly, for a L4 
(T e ff ~1820 K) object, the forsterite cloud base is found to be situated at 2.2 bar while 
Woitke & Helling (2004) found it to be at 2.7 bar for Ti0 2 cloud. 

Theoretical investigation (Tsuji 2004a) claims that the thickness of the dust clouds and 
hence the location of the cloud deck influences the spectral energy distribution of L and T 
dwarfs. Also, it is suggested that the vertical height of dust cloud may vary for a given T e ff 
and surface gravity (Knapp et al. 2004; Tsuji 2004a). The degree of polarization too is 
strongly dependent on the vertical height of the dust cloud and hence on the location of the 
cloud deck. In C03, the cloud deck is considered to be at one scale height above the base. 
The condensation curve P c j decreases exponentially with the decrease in the atmospheric 
temperature and the value of P c j become negligibly small at T = 1600fT. For different 
atmospheric models, this temperature is attained at different atmospheric pressure height. 
The position of the cloud deck for different spectral types is also presented in figure 2. From 
figure 2 we note that the thickness of the cloud decreases as one goes from L8 to L0. In other 
words, the thickness of the dust cloud decreases with the increase in effective temperature 
for a fixed value of surface gravity. For L dwarfs hotter than L2, the cloud is very thin, 
much less than one scale height. For a L8 object, the forsterite cloud deck is calculated to be 
situated at about 4 bar pressure height. For a similar atmospheric model, Woitke & Helling 
(2004) calculated the Ti0 2 cloud deck at 0.24 bar. For a L4 object the forsterite cloud deck 
in our model is at 0.5 bar pressure height while the TiC>2 cloud deck calculated by Woitke & 
Helling (2004) for a similar object is at 0.1 bar. Therefore the cloud thickness in our model 
is smaller than that of Woitke & Helling (2004). The number density of Ti0 2 grains for L4 
and L8 objects is presented graphically in Woitke & Helling (2004). 

At present, there is no convincing justification in favor of any specific form of the 
particle size distribution function. C03 considered a particle size distribution function that 
is consistent with the measurements of grain distribution attained in Earth's water clouds 
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while Ackerman & Marley (2001) and Saumon et al. (2000) considered a broad lognormal 
size distribution. In the present work, we adopt the size distribution function given by the 
latter authors which is expressed as : 



f(d) = ^-exp 



\n 2 (d/d ) 
In 2 a 



(20) 



where d is the diameter of the particles and d is the mean diameter, a is fixed at 1.3 that 
provides good fit to the red spectrum of L dwarfs. 



6. Results and discussion 

In a heterogeneously condensing cloud model with forsterite as the dominant constituent, 
the number density of grains, the cloud base and its vertical scale height are fixed. The 
remaining free parameters are the mean particle diameter d and the oblateness induced 
by rotation. As discussed by C03, the particle sizes significantly vary with the effective 
temperature, surface gravity and vertical height of the dust cloud. In the present work we 
have calculated the degree of polarization for a fixed surface gravity g = 10 5 cm s -1 . Figure 
4 shows how the degree of polarization is altered with different mean particle diameter 
for objects with different spectral type and hence for different effective temperature. In 
a multiple scattering scenario, it is important to consider different mean particle size at 
different pressure scale. However, in a single scattering, the photon scatters only with one 
particle and therefore, for a fixed effective temperature, a fixed mean particle size is sufficient. 
This means that at every atmospheric pressure level, we use the lognormal distribution with 
particle mean diameter d . 

We have included contributions to polarization by multipoles 1—2, 4 and 6 in eq. 14. 
However, it is found that at optical wavelengths and for the particle size needed to account 
the observed polarization, contribution from I — 2 is dominant over that of higher values for 
I. 

The projected rotational velocity for several L dwarfs has been determined from obser- 
vational data. In the absence of any knowledge on the inclination angle, the actual rotational 
velocity cannot be determined, yielding uncertainties in the oblateness. In the present work 
we have considered rotational velocity of 15, 25, 30, 40 and 50 km s _1 . Figure 5 and figure 6 
show how the degree of polarization significantly increases with the increase in rotational 
velocity. 

A comparison of the degree of polarization presented in figure 5 and figure 6 shows that 
if the mean particle diameter is increased, the degree of polarization decreases significantly 
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because with the increase in grain size, the grain number density reduces according to eq. 18. 

As mentioned in section 5, the location of the cloud base plays an important role in 
determining the degree of polarization. As one moves from LO to L8, the cloud base goes 
deeper in the atmosphere owing to the decrease in effective temperature. This makes the 
vertical size of the cloud larger (figure 2) and hence substantial increase in the total number 
of dust grains. In figure 3, we show how the degree of polarization is altered when the cloud 
base is changed but the particle number density remains unaltered. For an object of spectral 
type L4.0 (T c g ~ 1820 K) with surface gravity g = 10 5 cm s~ 2 , the cloud base should be 
situated at an atmospheric temperature about 1750 K. Keeping the particle number density 
the same as that calculated for L4 object, we calculated the degree of polarization by varying 
the location of the cloud base. The degree of polarization is found to increase as the location 
of the cloud base is pushed to the deeper in the atmosphere where the temperature increases 
from 1600 K to 1800 K. This implies that as one moves from L8 to L0, that is from cooler 
to hotter dwarfs, the degree of polarization decreases owing to the decrease in the vertical 
height of the dust cloud. 

As mentioned in section 1, Menard et al. (2002) for the first time, detected linear 
polarization from three L dwarfs. They use I bessel filter having the central wavelength at 
0.768 fim. Out of the three objects that show confirmed polarization, one object (DENIS-P 
J0255-4700) belongs to the spectral type L8 and has high rotational velocity (40 km s _1 ). 
This object shows a linear polarization of 0.167%. On the other hand, Zapatero Osorio et al 
(2004) have reported confirmed polarization from nine L dwarfs in I band (centered at 0.85 
\xm) and one L dwarf (2MASS J0036+1821) in I as well as in R band (centered at 0.641 \ira). 
The object 2MASS J0036+1821 is also observed by Menard et al. (2002) in I band but at 
different central wavelength. Hence, the degree of polarization at three different wavelength 
region is available for this object. 

Among the ten L dwarfs that are found to show confirmed polarization by Zapatero 
Osorio et al (2004), one (2MASS J2244+20) is reported to be significantly redder than 
those of other mid-L to late L dwarfs in near-infrared and infrared. Out of the remaining 
nine L dwarfs, the projected rotational velocity of only four objects is known. The object 
2MASS J2252-17 belong to the spectral type L 7.5. Although it's rotational velocity is not 
known, it shows degree of polarization as high as 0.62 %. Unless the rotational velocity of 
this object is very high (even higher than the fastest L dwarf Kelu 1), the large difference 
in the observed degree of polarization from this object and DENIS-P J0255-4700 cannot be 
explained as the mean particle size in the atmosphere of L dwarfs of same spectral type 
should not differ much. Further theoretical investigation on the distribution and location of 
condensates in the coolest L dwarfs are needed before modeling their degree of polarization. 
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Hence, in the present work we discuss linear polarization of L dwarfs of spectral type ranging 
from LO to L6. 

Figure 5 shows the degree of linear polarization at A = 0.850 /im for objects with 
different spectral types but fixed surface gravity 10 5 cm s -2 . It is found that the observed 
degree of polarization of several L dwarfs can be fitted if the mean diameter of grain is 
taken as 1.4 fim and with the polytropic index n=1.0. The figure shows that the change 
in polarization with the spectral type is not linear but overall the degree of polarization 
decreases as one moves from LI to L6. The degree of polarization decreases substantially 
for objects hotter than LI because condensation is not favored at such high temperature. 
The location of the cloud base shifts to deeper region of the atmosphere as the effective 
temperature decreases and hence should cause an increase in the polarization (figure 3). 
However, the total amount of condensing material is conserved. Therefore, as the scale height 
of the cloud layer becomes smaller, the particle number density becomes higher yielding 
into higher polarization. Beyond LI object, the temperature becomes too high to favor 
condensation and the polarization falls rapidly to zero. 

Zapatero Osorio et al (2004) detected the polarization in the Johnson R and I-band 
filters centered on 0.641 and 0.850 /im respectively with the passband of these filters as 0.158 
and 0.15 pm. Menard et al. (2002) detected the polarization in the Bessel I-band filter on 
0.768 fim with the passband as 0.138 \xm. We have calculated the degree of polarization at 
the central wavelengths. The change in degree of polarization over the spread of wavelength 
should by and large be absorbed in the error bars. 

Both 2MASS J1707+43 and 2MASS J1412+16 belong to the spectral type of L0.5. 
2MASS J1412+16 having projected rotational velocity 16.4 kms -1 shows degree of polar- 
ization 0.57±0.19 while 2MASS J1707+43 shows degree of polarization 0.23±0.06 but its 
rotational velocity is not known. Figure 5 shows that the observed data of these two objects 
can well be explained if the rotational velocity of 2MASS J1707+43 is much less than 15 
kms -1 . The model assumes the mean diameter of grain is 1.4 \im and the polytropic index 
n=1.0. The same model can explain the observed polarization of 2MASS J1507-16 (L5 with 
projected rotational velocity 27.2 kms" 1 ) if its effective temperature is slightly higher then 
that given by the spectral type - T e ff polynomial formula. We put it as L5.5 to show that if 
the effective temperature of this object corresponds to spectral type L5.5 in the polynomial 
formula then the observed polarization can be well fitted. The same model can explain the 
observed polarization from 2MASS J2158-15 (L4), 2MASS J0141+18 (L4.5) and 2MASS 
J0144-07 (L5) if their rotational velocities are about 45, 25, 25 kms -1 respectively. However, 
this model fails to fit the data from Kelu 1 (L2.5) unless its actual rotational velocity is 
between 25 to 30 kms -1 . 
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Figure 6 presents the degree of polarization with the same model but with larger mean 
particle diameter. It is found that the observed degree of polarization from 2MASS 1707+43 
(L0.5) can also be explained if its rotational velocity is 20 kms -1 but the atmosphere contains 
grains with mean diameter 3.0 \xm. However, the mean grain size of objects with the same 
spectral type should not differ much and therefore we predict the rotational velocity of this 
object is about 5-10 kms -1 with mean particle diameter 1.4 jim as presented in figure 5. 
Figure 6 shows that the observed polarization of Kelu 1 (L2.5) can be explained if the mean 
particle diameter is 3.0 jj,m. However, it is worth mentioning here that the actual rotational 
velocity of L dwarfs can not be determined from their projected rotational velocity unless the 
inclination angle is known. Figure 6 shows that the model with mean particle size 3.0 \xm and 
polytropic index n=1.5 ( non-relativistic completely degenerate polytropic distribution) can 
explain the observed polarization of 2MASS J0141+18 (L4.5) and 2MASS J0144-07 (L5.0) 
if their rotational velocities are the same of Kelu 1 (60 kms" 1 ). 

As mentioned earlier, the object 2MASS J0036+1821 (L3.5) is observed in three different 
wavelength regions and the degree of polarization is found to decrease substantially with the 
increase in wavelength. This trends strongly supports the presence of dust in the atmosphere 
of L dwarfs as it is very much unlikely that any other mechanisms such as the presence of 
magnetic field can explain this. Figure 7 presents the degree of polarization as a function 
of wavelength for L3.5 object with the polytropic index n=1.5 and the rotational velocity 
v=15 kms" 1 . It is obvious from the observed data at three different wavelengths that the 
atmosphere of 2MASS J0036+1821 should have grains of sub-micron size. We find the best 
fit of the three observed data with the mean particle diameter d =0.43 /an. Larger grain 
size would have made the polarization to peak at longer wavelength. 

The degree of polarization increases with the decrease in the polytropic index because 
with n = 1.0, the oblateness of the object is higher than that with n = 1.5 for the same 
rotational velocity and surface gravity. On the other hand Sengupta (2003) showed that 
the oblateness decreases with the increase in surface gravity yielding into less amount of 
polarization. 

As mentioned before, we have considered only forsterite in our models as it is very 
common in the atmosphere of L type dwarfs with solar metalicity. This is because of large 
abundances of Mg, Si and O. However, gehlenite, enstatite etc. should also be present in 
the atmosphere in fairly good abundance. The condensation curve presented in C03 shows 
that the base of gehlenite is situated much deeper in the atmosphere than that of forsterite 
yielding into a larger verticle size of the cloud. Further, iron with higher refractivity may 
undercircle most of the silicate clouds. Inclusion of all these will lead to substantial increase 
in degree of polarization and therefore, much smaller grain size may be needed in order to 
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explain the observed degree of polarization. 

Lastly, multiple scattering will lead to much less polarization and therefore in order to fit 
the observed polarization one has to either consider much higher oblateness of the objects or 
has to increase the grain number density substantially. An increase in grain number density 
needs smaller grain size which may contradict the present theoretical understanding on the 
nature and formation of dust grain. On the other hand, differential photometric observation 
of several L dwarfs could not detect any non-periodic variability of many objects that show 
high polarization. This indicates optically thin dust shell and therefore polarization by single 
dust scattering is quite reasonable. 

7. Conclusions 

We have investigated the optical linear polarization from L dwarfs of different spectral 
type by considering single dust scattering. Forsterite is considered to be the dominant species 
among the various condensates that could be present in the atmosphere of L dwarfs. The 
location of the cloud base and the cloud deck is determined from the condensation curve for 
forsterite and the atmospheric temperature-pressure profiles of different spectral types. The 
surface gravity is fixed at 10 5 cm s -2 and a wide range of rotational velocity is considered. 
It is found that the degree of linear polarization decreases from hotter to cooler L dwarfs. 
However, L dwarfs with effective temperature greater than 2200 K should not show detectable 
amount of polarization due to dust scattering as most of the dust would either evaporate 
from the atmosphere or condensation is not favored at such high temperature. It is found 
that the mean diameter of grains that is consistent with the observed polarization should not 
exceed a few micron although a small amount of very large grains at the base of the cloud 
for comparatively cooler L dwarfs may well be accommodated. However, the observational 
data of 2MASS J0036+1821 clearly indicates the presence of sub-micron size grain. Further 
polarimetric observation at the optical and at other wavelengths would provide convincing 
information on the properties and distribution of dust in the atmosphere of L dwarfs. 

We are thankful to the referee for several valuable comments and suggestions. 
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Fig. 1. — Effective temperature vs spectral type : solid line represents the calibration of T e ff 
using the sixth order polynomial fit as given in Golimowski et al. (2004) while the dashed 
line represents that by using the linear relationship given in Stephens et al. (2001). 
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Fig. 2. — Locations of forsterite cloud base (upper curve) and cloud deck (lower curve) for 
different spectral type L. 
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Fig. 3. — Degree of polarization for L4.0 object with the cloud base situated at different 
temperature in the atmosphere but with a fixed particle number density. The star indicates 
the degree of polarization when the cloud base is situated at 1750 K as determined from the 
intersection of the condensation curve and the T-P curve for L4.0 with g = 10 5 cm s~ 2 . 
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Fig. 4. — Degree of polarization as a function of mean grain diameter do for L dwarfs with 
different spectral type (L1-L6). 
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Fig. 5. — Degree of polarization for objects with different spectral type with the mean grain 
diameter rf = 1-4/xm and the polytropic index n=1.0. The five curves are for rotational 
velocities of V=15, 25, 30, 40 and 50 km s -1 , from bottom to top. The observed polarization 
of seven L dwarfs are plotted with their respective error bars. 
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Fig. 6. — Same as figure 4 but with the mean grain diameter rf = 3.0//m. The solid curves 
represent the degree of polarization with the polytropic index n=1.0 while the dashed curves 
represent that with n=1.5. The four pair of curves are for rotational velocities of V=20, 30, 
40 and 60 km s _1 from bottom to top. The observed polarization of four L dwarfs are also 
presented. 
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Fig. 7. — Degree of polarization as a function of wavelength with the mean grain diameter 
do = 0.43/im and polytropic index n=1.5. The observed polarization of the L dwarf 2MASS 
J0036+1821 at three different wavelengths are fitted with this model. The rotational velocity 
of this object is taken as V=15 km s _1 



